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ABSTRACT 

Due to its proximity, the mass of the supermassive black hole in the nucleus 
of Andromeda galaxy (M31), the most massive black hole in the Local Group 
of galaxies, has been measured by several methods involving the kinematics of 
a stellar disk that surrounds it. We report here the discovery of an eccentric 
Ha emitting disk around the black hole at the center of M31 and show how 
modeling this disk can provide an independent determination of the mass of the 
black hole. Our model implies a mass of 5.0^'q x 1O 7 M for the central black 
hole, consistent with the average of determinations by methods involving stellar 
dynamics, and compatible (at la level) with measurements obtained from the 
most detailed models of the stellar disk around the central black hole. This 
value is also consistent with the M — a relation. In order to make a comparison, 
we applied our simulation on the stellar kinematics in the nucleus of M31 and 
concluded that the parameters obtained for the stellar disk are not formally 
compatible with the parameters obtained for the Ha emitting disk. This result 
suggests that the stellar and the Ha emitting disks are intrinsically different from 
each other. A plausible explanation is that the Ha emission is associated with a 
gaseous disk. This hypothesis is supported by the detection of traces of weaker 
nebular lines in the nuclear region of M31. However, we cannot exclude the 
possibility that the Ha emission is, at least partially, generated by stars. 

Subject headings: galaxies: nuclei — galaxies: individual(M31) — galaxies: kine- 
matics and dynamics — techniques: spectroscopic 



Introduction 



All massive galaxies appear to host a supermassive bl ack hole (with 10 — 10 M c 



at their center ( IGultekin et al.l l2009t iMcConnell et al.ll201l[ ). Measuring the mass of cen- 
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tral black holes in galaxies is of great importance, as the discovery of a relationship be- 
tween mass and the veloci t y dispersion of the sta rs in the central bulge, the M — a relation 



f lFerrarese fc Merrittll2000l ; iGebhardt et al 



holes and their host galaxies (IGranato et al.l 120041 ) 



2000), reveals the possible co-evolution of black 



M31, the Andromeda galaxy, is an Sb galaxy at a distanc e of 77 8 kpc and its nucleus 
can be observed with excellent spatial resolutions. iLight et al.l ( 119741 ). using data obtained 
with the Stratoscope II, revealed an asymmetry in the nuclear region of M31, as the bright 
nucleus did not coincide with either the center of the bulge or the maximum of the stellar 
velocity dispersion. However, lLauer et al.l (119931 ). using observations from the Hubble Space 
Telescope (HST), showed that the galaxy possesses a double nucleus, the two components 
being called PI (the brightest one) and P2 (located, approximately, at the center of the 
bulge). These two components are separated by about 0'.'49. 



A model to explain the morphology of the nucleus of M31 was proposed by iTremaine 



( 119951 ) and states that PI and P2 are parts of an eccentric stellar disk around the black 
hole, with PI coinciding with the apocenter and t he black hole being located at P2. Several 
refinements to this model have been put forth (jSalow fc Statlerl l2004t iPeiris fc Tremaine 
20031 ); iBender et al.l ( 120051 ). using HST data, revealed that the black hole is actually located 
in a st ructure emb e dded in P2 called P3, which probably corresponds to a cluster of A- type 
stars. lLauer et al.l (120121 ). using also HST data, confirmed that P3 corresponds to a cluster 



of blue stars around the central black hole. 

The mass of the central black hole of M31 has already been measured by, at least, six dif- 
feren t techniques : (1) s tandard dynamical modeling ignoring asymmetries (IDressler fc Richstone 
19881 ; lKormendvlll988l); (2) the center of mass argument, which depends on the asymmetry 
of P1+P2 ( iKormendy fc Benderlll999l ); (3) dyn amical modeling; of the stellar nuclear disk 
taking into account the asymmetry of P1+P2 (IPeiris fc Tremaind 120031 ); (4) complete dy- 
namical modeling taki ng into account the a symmetries and the self-gravity of the nuclear 
stellar disk of P1+ P2 (Salow &: Statlerl 120041 ); (5) dynamical mod eling; of P3, which is inde- 



pendent of P1+P2 (IBender et al.ll2005l ): (6) iV-body simulations dBacon et al.l 120011 ). All of 
these methods involved stellar dynamics and resulted in values in the range 0.3 — 23 x 10 7 M Q 
for the mass of the central black hole in M31. 

In this Letter, we analyze a data cube of the nuclear region of M31, obtained with the 
Integral Field Unity (IFU) of the Gemini Multi-Object Spectrograph (GMOS) of the Gemini 
North telescope, and report the discovery of an eccentric Ha emitting disk around the central 
black hole. 
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2. Observations, reduction and data treatment 

The observations of M31 were made on 2009 September 21. We used the IFU of the 
GMOS of the Gemini North telescope, in the one-slit mode, in order to obtain data cubes, 
with two spatial dimensions and one spectral dimension. The science field of view (FOV) 
has 5" x 3'.'5, while the sky FOV (observed simultaneously at a distance of V from the science 
FOV) has 5" x 1".75. Three 10 minute exposure of the nuclear region of M31 were made, 
with the grating B600-G5307, in a central wavelength of 6000A. The final spectra had a 
coverage of 4550 — 7415A and a resolution of R ~ 2900. The estimated seeing for the night 
of observation was 0'.'55. 

Standard calibration images were obtained during the observations. The data reduction 
was made in IRAF environment. At the end of the process, three data cubes were obtained, 
with spaxels of 0'.'05 x 0'.'05. No sky subtraction was applied because the sky FOV (still inside 
the disk of M31) was contaminated with stellar emission from the galaxy. 

After the data reduction, we performed a procedure of data treatment. First, a cor- 
rection of the differential atmospheric refraction was applied to all data cubes, using an 
algorithm developed by our group. In order to combine the three corrected data cubes into 
one, a median of these dat a cubes was calculated. After that, a Butterworth spatial filtering 



(IGonzalez fc Woods! 120021 ). with order n = 2, was applied to all the images of the resulting 
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using a synthetic Gaussian point-spread function (PSF). The PSF of the final data cube has 
FWHM ~ 0'.'45. 

Figure [1] shows an image of the final data cube of M31 (obtained after the data treat- 
ment) collapsed along the spectral axis and an average spectrum of this data cube. The 
brightest component of the nucleus, PI, can be easily detected; however, the fainter compo- 
nents, P2 and P3, cannot be seen, due to the spatial resolution and to the lack of spectral 
sensitivity in the blue (below 5000 A). A spectrum of PI, extracted from a circular area with 
a radius of 0'.'15, is also shown in Figured! The average signal-to-noise ratio (S/N), between 
5610A and 5670A, of the spectra of the data cube analyzed here is close to 50. 



3. Data analysis and results 



After the data treatment, a spectral synthesis was applied to the spectrum of each 
spaxel of the resulting data cube of M31. This procedure was performed with the Starlight 
software (jCid Fernandes et al.l 120051 ) . which fits the stellar spectrum of a given object with 
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a combination of template stellar spectra from a pre-established base. In this work, we used 
the base of stellar spectra MILES (Medium resolution INT Library of Empirical Spectra; 
Sanchez-Blazquez et al. 2006). The spectral synthesis resulted in a synthetic stellar spectrum 
for each spaxel. These synthetic spectra were then subtracted from the observed ones, leaving 
a data cube with emission lines only. The non subtraction of the sky field during the data 
reduction had no observable effect in the results obtained with the spectral synthesis. 

In this residual data cube, we have detected a previously unreported weak extended 
Ha emission. The luminosity of this Ha emission line, within a radius of 0".7 from the cen- 
tral black hole, is Lf/ a (emission) ~ (8.7 ± 1.0) x 10 2 L Q ; the equivalent width of this line is 
WHa(emission) m — 0.18±0.03A, while that of the absorption component is WHa(stellar) « 
2.03±0.10A, indicating that the extended Ha emission is much weaker than the stellar com- 
ponent and, therefore, hard to detect. Figure [2] shows the spectral fit of the Ha absorption 
line in one spectrum of the data cube of M31 and also the fit residuals. It is possible to 
detect traces of an Ha emission line immersed in the absorption component, however, only 
the subtraction of the spectral fit allowed a clear visualization of this very weak emission 
line. Figure |2] also reveals the presence of H/3 and [O III] A5007 emission lines. The spectral 
resolution of the base of stellar spectra MILES (FWHM = 2.3A) is very similar to our spec- 
tral resolution (FWHM = 2.2A) and, because of that, the spectral fit is quite precise and 
allows an accurate subtraction of the stellar emission (Figure [2]). The spectra of the data 
cube of M31 were resampled to AA = I A (the same spectral sampling of MILES spectra), 
before the spectral synthesis was applied. 

In Figure |2j we can see a map of the integrated flux of the Ha emission line. Only 
the central region of the FOV of the data cube was taken into account here, because the 
Ha emission was too faint farther away. We can see two bright areas in the vicinity of PI 
and P2, respectively. This reveals a certain similarity between the extended Ha emission 
and the double nucleus of M31 and suggests the existence of a possible Ha emitting disk. 
An elongated emitting region from SE toward NW can also be seen; however, we do not 
know the origin of this peculiarity in the morphology of the extended Ha emission. The 
area wit h masked values in th e map corresponds to an emitting region, observed in previous 



studies f ldel Burgo et al.ll2000l ). that is not associated to the extended Ha emission detected 
here. This emitting region can be identified by observing its intense, compact, and narrow 
[O III]AA4959,5007 emission. 

In order to obtain a velocity map for Ha, we fitted a Gaussian function to the Ha 
emission line in each spectrum of the data cube. Figure [2] shows an example of the Gaussian 
fit applied to the Ha emission line detected in one spectrum of the data cube, after the 
subtraction of the spectral fit. We can see that, despite the irregularities of the emission 
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line, the Gaussian fit provides (with a considerable precision) the wavelength corresponding 
to the peak of the emission line and, therefore, the radial velocity. On the other hand, the 
irregularities observed in the Ha emission lines, after the subtraction of the spectral fit, made 
impossible to determine reliable values for the velocity dispersion. Therefore, such values 
were not taken into account in this work. 

The velocity map obtained for the Ha emission line is shown in Figure |3] and it clearly 
indicates the presence of an Ha emitting disk around the central black hole. The position 
angle of the line of nodes of this emitting disk is P.A. e = 59?3 ± 2?7, which is very close to 



the y alue measured for the stellar disk in previous studies (P.A.j = 56?4±0?2) (IBacon et al. 



20011 ). The contours in Figure [3] show that the ascending node of the Ha emitting disk is 
very close to PI, while P2 and P3 are located in an area with low velocities. This behavior 
is very similar to what is observed in the stellar disk. We extracted a velocity curve along 
the line of nodes (Figure [3]) and the modulus of the maximum and minimum velocities are 
considerably different, indicating that any kinematical model should be eccentric. 

We tried to reproduce the velocity curve and the velocity map of Ha using a model of 
a simple eccentric disk around the supermassive black hole. Only a region within a radius 
of 0".7 from the black hole was modeled because, as mentioned before, the disk was too faint 
farther away. We admitted that the stellar mass inside the radius of the modeled area was 
small compared to the mass of the black hole, so no stellar mass was taken into account 
in the model (see more details below). The emitting disk was simulated by superposing 33 
concentric Keplerian orbits with the following free parameters: the argument of the peri- 
center u, the inclination of the disk i, the eccentricity of the disk e and the mass of the 
central black hole M,. We measured the value of the longitude of the ascending node in the 
velocity map (Q = 128?3), so this parameter was fixed in our model. In each simulation, 
after all the orbits were superposed, the model was convolved with the estimated PSF, in 
order to simulate the effect of the Earth's atmosphere. The free parameters were varied and 
the simulations were made repeatedly, in order to minimize the value of the y 2 between the 
observed velocity map and the simulated one: 

2 1 Wij ■ Iij ■ (vij (observed) — Vij(simulated)) 2 



7 EE 



X - , / j ^ 2 

1 j=l a h 

where N x is the number of spaxels along the horizontal axis, N y is the number of spaxels 
along the vertical axis, is the uncertainty of the velocity of the spaxel (i, j), Vij(observed) 
is the velocity of the spaxel (i,j) of the observed velocity map, i>y (simulated) is the velocity 
of the spaxel (i,j) of the simulated velocity map, is weight equal to 1 for areas near the 
line of nodes and equal to for farther areas, ly is the value of the spaxel (i,j) of the map 
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of integrated fluxes of Ha (Figure [2]), and / is the sum of the values of all spaxels of the map 
of integrated fluxes of Ha (Figure [2]). 

In the previous formula, the weight Wy is represented by a step function, which is equal 
to 1 for spaxels closer than 0'.'25 to the line of nodes and equal to for farther spaxels. We 
used the values of the integrated flux of the Ha emission line (i^) in the calculation of the 
X 2 in order to give more weight to the spaxels with higher S/N. 

The best simulated velocity map and the corresponding best simulated velocity curve 
are also shown in Figure [31 The parameters of the best model obtained, with \ 2 = 1.67, 
are shown in Table [TJ The uncertainties (la) of the parameters were estimated by plotting 
histograms (probability distributions) of each one of the parameters of the simulation, con- 
sidering, only cases in which X 2 ~Xmin < 1- After that, we fitted a Gaussian function on each 
histogram and obtained the sq uare deviation, which was taken as the uncertainty of each 



parameter ( jde Rosa et al.ll201ll ). We can see that, despite some irregularities of the observed 
velocity map, our model of a simple eccentric disk reproduced the observed kinematical be- 
havior considerably well. In Figure |3] however, we can see a clear discrepancy between the 
observed velocity map and the simulated one in an area near the ascending node. This is the 
same emitting region that was masked in the map of the integrated flux of the Ha emission 
line (Figure which is not associated to the eccentric disk detected here. A by-product of 
this modeling was an independent determination of the position of the black hole. We found 
that it is at a distance of 0'.'050 ± 0'.'025 from P2, which is compatible (at la le vel) with the 



posit ion of P3 (at a distance of 0'.'076 from P2) from previous determinations (IBacon et al. 



20011)- 



As mentioned before, in this model, we have not taken into account the effect of the mass 
of the stars. In order to evaluate this assumption, we used an HST image of the nucleus of 
M31 in V band, obtained with WFPC2, to estimate the stellar mass within the radius of the 
simulated area. First, we decomposed this image into an asymmetric component (containing 
the stellar disk around the black hole) and a symmetric one (containing the central part of the 
stellar bulge). After t hat, we de-projec ted the symmetric component using a multi-Gaussian 



expansion approach ( ICappellaril 120021 ) . Finally, the stella r masses of the two components 



were calculated using a mass-to-light ratio of M/L v = 5.7 (IBender et al.l 120051 ) . The results 
are M sym = 8.7 x 1O 5 M for the bulge component and M asym = 8.7 x 1O 6 M for the disk. 
Together, these two components represent 19% of the M. we obtained. That means that 
this value of M, resulting from our model can be, at most, 19% too high, due to the neglect 
of the stellar mass. Introducing this in the previous uncertainty for the black hole mass, we 
obtain a final value of Mbh_Ha = 5.0l^o x 1O 7 M . Since the spectral synthesis performed on 
the data cube provided details about the stellar populations detected, we could, in principle, 
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have used a mass-to-light ratio derived from these results to estimate the stellar mass in the 
nuclear regio n of M31. However, we decided to use a mass-to-light ratio derived from HST 



observations (IBender et al.ll2005l ) because, due to our spatial resolution and lack of spectral 
sensitivity in the blue (below 5000 A), we did not observe certain young stellar populations 
in the nucleus of M31 that were clearly detected in HST observations. 

We performed an entirely analogous modeling on a stellar velocity map obtained from 
the data cube of M31. Only a region within a radius of 0'.'8 from the black hole was modeled 
because the effect of the self-gravity of the stellar populations could be considerable at farther 
areas. The parameters of the best model obtained are shown in Table [1] and the final mass 
obtained for the central black hole (including the uncertainty introduced by the neglect of 
the stellar mass in the modeling) is M bh _ ste u ar = 4.515!? x 1O 7 M . We will show the details 
about this topic in R. B. Menezes et al. (2013, in preparation), since the stellar kinematics 
analysis is not the purpose of this Letter. 



Discussion and conclusions 



The mass of the central black hole in M31 we obtained by modeling the kinematics 
of the Ha emission as a simple eccentric disk (Mbh_Ha = 5.0jZ°'o x 1O 7 M ) is within the 
ra nge of the v a lues f ound in previous dete rminations using stellar kine matics (10 6 * 5 — 1O 8 M 



in 



Kormendvl (119881): 3 



Dressier & Richstone ( 


1988); 


o in 


Peiris & Tremaine 


2003) 



Bender et all (120051 ): 7.0±3;5 x 1O 7 M in Bacon et al.l (1200 ll )) and it is compatible (at la 



-1.2 



x 1O 7 M in 



level) with the measurement (5.62 ± 0.66 x 10 7 M^) obtaine d with the most detailed models 



of the stellar disk elaborated so far (jSalow fc Statlerl 12004 ). This result is also compatible 



(at la level) with the value of the mass of the central black hole obtained by modeling the 
stellar kinematics as a simple eccentric disk (Mt,h_ s teUar — 4.51?'? x 1O 7 M ). Most of the 
recent models used to reproduce the stellar disk around the black hole in M31 did not use a 
single inclination and eccentricity for the entire disk. However, the values we obtained for % 
and e in our best models are in the range of the ones used in several stud ies to model the 
stellar disk (0 — 0.7 for the eccentricity and 4 1° — 77° for the inclination) (jSalow fc Statler 
2004 : iTremaind Il995l ; iPeiris fc Tremaind 120031 ) . 



The values of e and i obtained for the stellar and the Ha emitting disks with our models 
are compatible at la level. This reveals a certain similarity between these two models. 
However, the values of u obtained are not compatible at all. This suggests that the stellar 
and the Ha emitting disks are intrinsically different from each other. According to what 
was mentioned before, we have identified a weaker H/3 emission line and also traces of even 
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weaker nebular lines, like [N II] A6583, [S II] AA6716, 6731, and [O III] AA4959, 5007 (Figure^, 
along the FOV of the residual data cube of M31. This, together with the intrinsic difference 
between the stellar and the Ha emitting disks, suggests that the Ha emission is associated 
with a gaseous disk. However, we cannot exclude the possibility that this emission is, at 
least in part, originated by stars. 

The discovery of an Ha emitting disk in the nucleus of M31, reported here, allows an 
independent measurement of the mass of the central black hole in M31, and, therefore, has 
a considerable importance for the studies of this object. 

This work is based on observations obtained at the Gemini Observatory. We thank 
FAPESP for support under grants 2008/11087-1 (R.B.M.) and 2008/06988-0 (T.V.R) and 
also an anonymous referee for valuable comments about this Letter. 

Facilities: Gemini: Gillett(GMOS). 
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Fig. 1. — Left above: image of the final data cube of M31 collapsed along the spectral axis. 
Right above: average spectrum of the final data cube of M31. Bottom: spectrum extracted 
from a circular area, with a radius of 0'.'15, centered on PI. 
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Fig. 2. — a) Ha absorption line of one spectrum of the data cube of M31 (+0'.'l;— 0'.'3), with 
the spectral fit (in red) and the fit residuals (in green), b) Same fit residuals shown in a), 
with a Gaussian fit (in red), c) [O III] A5007 emission line of the same spectrum shown in 
a), with the spectral fit (in red) and the fit residuals (in green), d) H/5 absorption line of the 
same spectrum shown in a), with the spectral fit (in red) and the fit residuals (in green), e) 
Map of the integrated flux of the Ha emission line. The positions of PI, P2, and P3 and the 
contours are also shown. 
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Fig. 3. — Left above: observed velocity map of the Ha line. Right above: best simulated 
velocity map of the Ha line, with the position of BH (the black hole in our model) indicated. 
The contours of the integrated flux of the Ha line (Figure E]) and the positions of PI, P2, 
and P3 are shown in the observed and in the best simulated velocity maps of Ha. Bottom: 
velocity curve of the Ha emission line, measured along the line of nodes. The best simulated 
velocity curve (in red) is superposed on the measurements (points with error bars). 
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Table 1: Parameters of the Best Model Obtained for the Ha Emitting Disk and for the Stellar 
Disk. 

I 1 ■ , : i n< ' - ■ Ha Emitting Disk Stellar Disk 

M. (5.0 ± 0.8) x 1O 7 M (4.5 ± 0.9) x 1O 7 M 

e 0.35 ±0.07 0.30 ±0.12 

u 115° ±12° 250° ±17° 

i 45° ±4° 55° ±10° 



